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^ ' A model with interlayer pairing is proposed to explain the sinusoidal modula- 

^ ■ tion of the resonant neutron scattering in high temperature superconductors. 

'^ . It is found that the interlayer pairing has s-wave symmetry in the Cu02 plane 

j^ , and has comparable magnitude with the d-wave intralayer pairing. It is also 
found that the interlayer pairing mainly affects momentum close to the hot 

__i - spots on the Fermi surface while its effect on the gap nodes is negligible. It 

Ch ' is pointed out that these characteristics of the interlayer pairing can be un- 
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derstood in a model in which the superconducting pairing originates from the 
exchange of the antiferromagnetic spin fluctuation. 



The discovery of the resonant neutron scattering is one of the most important progress 
in the high-T^ field in recent years |l|-§|. This resonance, also called (tt, tt) resonance, 
has many interesting characteristics and has attracted much theoretical attention 



The c-axis modulation is an interesting and intriguing problem in the explanation of the 
(tt, tt) resonance. In experiments, only the odd channel magnetic response [qc = vr) has been 
observed in the Cu02 bilayer systems (the sinusoidal modulation), while in most single Cu02 
plane-based RPA theories, the interlayer exchange is far too small to distinguish the even 



from the odd channel |rT|- [T3| , p!6| J]l 

Physically, the momentum dependence of the magnetic response is closely related to the 
internal structure of the superconducting Cooper pair. For point-like s-wave pairing, the 
magnetic response is generally momentum independent and is suppressed due to the singlet 
nature of the pair. While for d-wave pairing which takes place between nearest neighboring 
sites on a plane, the magnetic response is strongly momentum dependent. In fact, here we 
can take the Cooper pair roughly as the coherent superposition of two antiferromagnetic spin 
configuration. Thus, if we look at the pair with a momentum transfer of the order (tt, tt), we 
will find an enhanced magnetic response. Similarly, if there exists interlayer pairing between 
the two Cu02 planes in the Cu02 bilayer, the odd channel magnetic response {qc = vr) will 
be enhanced while the even channel response {qc = 0) will be suppressed. Although such 
pairing is obviously favored by the interlayer exchange coupling, it is totally neglected in 
the single Cu02 plane-based theories. 

In this letter, we find the sinusoidal modulation of the (tt, vr) resonance can be naturally 
explained with the inclusion of the interlayer pairing. We find that the interlayer pairing 
has s-wave symmetry in the Cu02 plane and has comparable magnitude with the intralayer 
d-wave pairing. We find that the interlayer pairing mainly affects momentum close to the 
hot spots on the Fermi surface and has negligible effect on the gap nodes. 

To model the Cu02 bilayer with interlayer pairing, we take the following mean field 
Hamiltonian [|T^ 

U - \^ C >)t>) I \^(A >)t>)t Mh r-\ m\^\a' |'^(l)t J2)t , J2)t (l)tN , 7, ^ 1 /I ^ 

^MF - 2^ t.kCka Cka + l^K^kCk^ C_fc^ + ll.C.) + }_^ [A^(Cfc^ C_f,^ + C^^ C_fc^j + tl.C.\ [1) 

k,n,<T k,n k 

in which n=l,2 is the layer index. C,k is the dispersion in the Cu02 plane. Here we use 
the dispersion derived from fitting the ARPES result in Bi2Sr2CaCu208 [ll8| , pO| , pT| , ^k = 



— t {cos kx + cos ky)—t' coskx cos ky —t"{cos2kx + cos2ky)—t"'{cos2kx cos ky + cos k^ cos2ky) — 
t""cos2k^cos2ky - fx, t = 0.2975 eV, t' = -0.1636 eV, t" = 0.02595 eV, t'" = 0.05585 eV, 
t"" = —0.0510 eV, fi is the chemical potential. Note that we have neglected the inter- 
layer hopping term in the dispersion since no band splitting is observed in experiment. 
Afc = ^(cos/Cx' ~ cos ky) is the intralayer d-wave pairing function. A^ is the interlayer 
pairing function. As will be shown later, A^ has strong momentum dependence in the Cu02 
plane. However, such momentum dependence is not essential for the discussion of the (vr, vr) 
resonance since the low energy magnetic response at g = (tt, tt) is determined mainly by the 
electronic transition between the hot spots on the Fermi surface (see FIG. 1). For the sake of 
simplicity, we will take A^ as momentum independent for the moment. The relative phase 
between A^ and A^ is another important issue. In the absence of the interlayer hopping, A^. 
can be either real or purely imaginary to meet the requirement of the time reversal symme- 
try. However, a real A^^ will lead to different energy gap at momentum [k^, ky) and {ky, k^), 



since A^ and A^ have different momentum dependence in the Cu02 plane. Therefore Aj, 
must be purely imaginary. 

To discuss the c-axis modulation of the Cu02 bilayer system, it is convenient to use the 
bonding band and the anti-bonding band representation 1^,11,16] 
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Here b and a represent the bonding and the anti-bonding band respectively. In this repre- 
sentation, the mean field Hamiltonian reads, 

Hmf = E Ufci'i + E(Ai"^4f cit + h.c.) (3) 

k,a,a k,a 

a;=a,b is the band index, A^"^ = A^ + /(a)A^, where 

„/ X I +1, for a = b 
^^") = \-l, fora = a 

In the bonding and the anti-bonding band representation, the even and the odd channel mag- 
netic response come from the intraband and the interband electronic transition respectively 



Xo"^'^^ (g, uj) = Xo''^ (g, ^) + Xo^^ {q, ^) 
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in which the mean field susceptibility Xo {Q)^)i'Xs) {Qi^)iXo {q^^) ^iid x\) {Qi^) ^i-re 
given by ( for simphcity we discuss the zero temperature case) ||lTl , p2| 



xV\q:^) = 7E(1 T^HTTJT) -)( . ^m . 



4. E^E^q "u + E^ + E^l + ^5 u - E^ - e\:1 + ^5' 

(5) 



here i, j=a, b, E]^' = \ {^k) + A^*^ is the quasiparticle energy. 

To see the effect of the interlayer pairing on the momentum dependence of the magnetic 

response, let us look at the BCS coherence factor, (1 Tir~&r~^^)' which contains the 

k fe + g 

information about the internal structure of the Cooper pair. In the absence of the interlayer 
pairing, the even and the odd channel have the same coherence factor (1 — "^'' ^ *" ''''"'' ) 

and both channels are fully enhanced at momentum transfer q = (vr, vr) since AkAk+q < 0. 
That is, the system dose not distinguish the even and the odd channel at the mean field 
level. As we have mentioned at the beginning of this paper, the RPA correction from the 
interlayer exchange is far too small to make this mean field result agree with the observed 
large difference between the even and the odd channel. 

In the presence of the interlayer pairing, the even and the odd channel behave differently. 

For the odd channel, since A^"^A^* = A^, A^!^* = — ( A| + A^ j (using the properties 



^k+q = — Afe, A).^ = A^), the coherence factor can still reach its maximum value 2 on the 
Fermi surface. That is, the odd channel magnetic response is still fully enhanced. While 



for the even channel, since A^ A^",:* = — ( A| — 
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A^ A^'^ (the cross term of A^ and A), vanishes upon summing over k since A^ and A). 
have different symmetry), the coherence factor can take any value ranging from (totally 



suppressed) to 2 (fully enhanced) on the Fermi surface depending on the ratio 



As a 



Ao 

result, the even channel magnetic response is suppressed with the increase of the interlayer 

pairing. FIG. 2. shows the calculated susceptibility for ^-^ = l(the magnitude of the 
interlayer pairing will be discussed later). We see the interlayer pairing suppresses the even 
channel magnetic response very effectively. 

So far, we have discussed the bare susceptibility xo- To obtain the fully renormalized sus- 
ceptibility, we still have to include the RPA correction from the antiferromagnetic exchange 
coupling. In the presence of the interlayer exchange coupling, the even and odd channel 
magnetic responsee are renormalized differently |[10|JTT|JT6|| , 



{even) , 
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in which Jq = J{cosqx + cosg^) is the intralayer exchange, and Jp is the interlayer exchange 
coupling. Experimentally, J ~0.15 eV, Jp/J ~ 0.1 P3|. FIG. 3 shows the renormalized 



susceptibility. After the RPA correction, a sharp resonance appears well below the gap edge 
in the odd channel [|l^]. While in the even channel, there is only a small peak very close to 
the gap edge. This result can be understood by examining the resonance condition for both 
channels. As can be seen from FIG. 2, the resonance condition is fulfilled well below the 
gap edge in the odd channel. While in the even channel, this condition is only fulfilled very 
close to the gap edge because of the reduced magnitude of the bare susceptibility in the even 
channel( Jp alone is too small to produce the observed even-odd difference). Here we find 
the magnetic response starts at different energies in the odd and even channel. This agrees 
very well with experimental observations [^ . According to our theory, the resonance energy 
of the odd channel is unrelated to the superconducting gap while the energy threshold for 
even channel magnetic response is about twice of the superconducting gap. 

In the foregoing discussion, we have neglected the momentum dependence of the inter- 
layer pairing. This is reasonable for the discussion of the (vr, vr) resonance which mainly 
concerns the hot spots on the Fermi surface. However, a momentum independent inter- 
layer pairing is inconsistent with the experimental observation of the gap nodes along the 
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To be consistent 



(0, 0) — (tt, tt) direction since the total energy gap equals y A| 

with the existence of the gap nodes. A;, must be negligibly small along the node direction. 
Thus, the interlayer pairing must be strongly momentum dependent in the Cu02 plane. 
Here, a closely related problem is the magnitude of the interlayer pairing. In our calcula- 

K 

tion, we have assumed -^— = 1. This may seem arbitrary at first sight. However, if we 
assume that the superconducting pairing originates from the exchange of the antiferromag- 



netic spin fluctuation |2^j2^, especially, the (vr, vr) resonance |^, then both the magnitude 



and the momentum dependence of the interlayer pairing can be easily understood. Since 
the resonance occurs only in the odd channel, the spin fluctuation mediating the intralayer 
and the interlayer pairing have the same propagator except for an overall sign( note that 

(even) _ (intralayer) , (interlayer) y(odd) _ (intralayer) _ (interlay er) \ Hgnce it is QUite 

reasonable that the intralayer and the interlayer pairing have comparable magnitudes (but 
different symmetry). At the same time, since the (vr, vr) resonance is sharply peaked at (tt, tt) 
IQJ^, only momentum close to the hot spots is significantly affected by the interlayer pairing. 
Hence the interlayer pairing must be strongly momentum dependent. 

Interestingly, this pairing mechanism also naturally explains the different symmetry of 
the intralayer and the interlayer pairing. This difference comes from the overall sign change 
between the intralayer and the interlayer spin fluctuation propagator. Since the exchange 
of intralayer antiferromagnetic spin fluctuation favors d-wave intralayer pairing p5|,^, or, 
^fc+(7r,7r) = —^k, the interlayer pairing mediated by the interlayer spin fluctuation must 
have s-wave symmetry, or, Af^^,. = A^^. When A^ is transformed into the real space, 
we will see the interlayer pairing exists only between sites of the same magnetic sublattice 
on the two Cu02 plane. s According to our discussion concerning the relation between the 
momentum dependence of the magnetic response and the internal structure of the Cooper 
pair, such interlayer pairing will enhance the odd channel response and suppress the even 
channel response, as we have observed. 

In conclusion, we find the sinusoidal modulation of the (vr, vr) resonance observed in 
experiments can be explained with the inclusion of the interlayer pairing in the theory. We 
find the interlayer pairing has s-wave symmetry in the Cu02 plane and has comparable 
magnitude with the d-wave intralayer pairing. We also find the interlayer pairing has strong 
momentum dependence and mainly affects momentum close to the hot spots on the Fermi 
surface. We find these characteristics of the interlayer pairing can be understood in a model 
in which the superconducting pairing comes from the exchange of the antiferromagnetic spin 
fluctuation, especially, the (vr, vr) resonance. 

The author would like to thank T.K.Lee, X.G.Wen, S.C.Zhang and F.C.Zhang for their 
valuable comments and Z. Tesanovic for pointing out Ref. 19 to him.. 
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FIGURES 

FIG. 1. The Fermi surface and the important momentum in our discussion. VHS denotes 
the Van Hove singularity. 

FIG. 2. The imaginary(a) and the real part(b) of the bare susceptibility, for 
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35 meV. The Fermi surface is 34 meV above the VHS. Note both the maximal 



gap and the chemical potential are the same as that used in Ref. [18]. The intersections of 
the straight lines and the curves in (b) give the resonance energies in the odd and the even 



channel, for J=150 meV, Jp=0.1 J. 



FIG. 3. The susceptibility after the RPA correction, for J=150 meV, Jp=0.1 J. 
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